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Considering the density-density correlation function of a concentrated polymer 
solution with multicomponent solvent we find a phase transition due to the hetero- 
geneity of excluded volume constant. This new phase transition implies a strong 
enhancement of the scattered light intensity in the critical region, which can explain 
a recent experiment showing strong light scattering from a ternary polymer system 
consisting of polyethylene oxide (PEO) dissolved in nitroethane and Smethyl-pentane. 
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Physical properties of polymer solutions have attracted much attention in recent years 
Polymer solutions with multicomponent solvent are of particular current interest 
because some polymers require mixed solvents to reach their G-conditions |^ and such 
polymer solutions have fluctuating excluded volume constants as in the case of random 
heteropolymers in a single solvent Therefore, the investigation of polymer solutions with 
multicomponent solvent will help to understand the properties of random heteropolymers, 
including many biopolymers such as RNA, DNA, protein, etc, which play important roles 
in functions of biological systems. 

Light scattering is one of the main methods for the investigation of polymer solutions 
1^. Particularly, it allows the determination of the molecular weight of polymers in dilute 
solutions 0. Many papers have studied the light scattering of polymer solutions with mul- 
ticomponent solvents [0. However, most of these papers considered the scattering from 
dilute polymer solutions. In the present paper we consider the concentrated polymer solu- 
tions with multicomponent solvent. In this case a strong overlapping between the polymer 
chains exists and the excluded- volume interaction between the monomers is screened 
Although the solution is concentrated it still occupies a very small volume fraction of the 
system (semidilute regime). Such a situation was realized in a recent experiment by To, Kim 
and Choi (TKC) P| who studied light scattering from a ternary polymer system consisting 
of polyethylene oxide (PEO) dissolved in a binary solvent of nitroethane (NE) and Smethyl- 
pentane (MP) (PEO/NE/MP) and found an abnormal strong light scattering, which does 
not have any theoretical explanation. Using the concept of a fluctuating excluded-volume 
constant to calculate the density-density correlation function of a polymer solution with a 
binary solvent, we find a critical behavior for such a system when the fraction of one com- 
ponent of the solvent approaches a critical value. This new phase transition implies a strong 
enhancement of the scattered light intensity in the critical region, which can explain TKC's 
abnormal strong light scattering in PEO/NE/MP 0. 

Consider a solvent the components of which have the same refractive index. This means 
that the solvent will not scatter light by itself. However suppose that the components 
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differ from each other by their abihty to dissolve a given homopolymer. Now consider 
a solution of concentrated homopolymer s in such a multicomponent solvent. Since the 
excluded-volume interactions between the monomers depend not only on themselves but also 
on their environments, they will fluctuate from point to point because of the composition 
fluctuations of the solvent. In analogy with the random heteropolymer case |TD[, the energy 



of a polymer solution with a fluctuating excluded-volume constant, in terms of monomer 
concentration f^, can be written in Gaussian approximation in the form 



^ = i / dfc{r) 
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c{r). (1) 



Here 

c(f)=5:5(f-/?0 (2) 

i 

is the microscopic concentration of monomers, c =< c{r) > is the average concentration of 
monomers, a is their average size and v is the average value of the excluded volume constant 
with fi(r) being its fluctuating part. The summation in Eq. runs over all chains and 
monomers and Ri is the coordinate of the i-th monomer. As was mentioned above, the 
excluded-volume constant fluctuations are caused by the composition fluctuations of the 
solvent and by the fact that the excluded-volume constant has different values in different 
components of the solvent. Therefore we assume that the fluctuating part f i(r) is a Gaussian- 
distributed random function whose correlation function is given by 

< Vi{r)vi{r') >= B{f- r'), < Vi >= 0. (3) 

The Fourier transform of the correlation function B[r — r') is given by 

Here w is the variance (dispersion) of the excluded-volume constant and is the correlation 
length of the composition fluctuations. Note that far away from the critical temperature 
of the solvent mixture is very small (less than 2nm for NE/MP mixture |]rT|) and fluc- 



tuations can be assumed to be almost uncorrelated. When the solvent consists of only two 



components, as in the experiment 0, the average and the variance of the excluded volume 
constant can be estimated as follows 

V = <!>Vp + {l-<t>)vg, = <l>{l-<^){vg-Vpf. (5) 

Here Vp and Vg are the excluded volume constants in poor (MP) and in good (NE) compo- 
nents, respectively; $ and 1 — $ are the fractions of MP and NE in the solvent, respectively. 

It is well known that the intensity of scattered light is proportional to the Fourier trans- 
form of the density- density correlation function of the polymer solution f^, which is deter- 
mined as follows 

G{f,r') = -\< c{f)c{r') > -c^] . (6) 

C L J 

Here < .... > means a thermodynamic average using the energy of Eq.(|l]). Thus 

cG(f, ?) = aJY[ dR,[c{^c{r') - c^]exp , (7) 

where A is a normalization constant and the integration goes over all monomer coordinates. 
It is convenient to go to integration over the monomer concentration c{r) instead of Ri in 
(0) 1^. Doing this and calculating the Gaussian integral over c(r*), we find that the correla- 
tion function G{f,r'), representing vc^'^G{f,r') above with ^ being a{l2vc)^^^'^, satisfies the 
equation ||10[ 
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G{f,?) = 5{r-?). (8) 



The first two terms in describe the polymer solution in a homogeneous solvent and the 
third one describes the heterogeneity of the excluded- volume constant. 
It is easy to obtain the bare [vi = 0) correlation function from Eq.(|]) 

Go{<i) = (9) 

which means that ^ is the correlation length of the polymer in a homogeneous solution. For 
averaging over the realizations of the random field f i , we use the impurity diagram technique 



1^ . Expanding Eq.(P) on Vi and averaging in each order over Vi, one can represent the 



averaged correlation function in a Dyson form 

where the self-energy S(g) is determined by the following irreducible diagrams 



(11) 



The heavy line denotes the bare correlation function of Eq.(^ and the dashed line is the 
Fourier transform of the random field Vi correlation function of Eq.(|^), C,^'^B{k)/v'^ . In the 
leading order we find from Eq. (p!l|) 



m = ^-^jj^amGo{\k-q\). (12) 



It follows from Eqs. (p!oD -(p!2D that fluctuations of the excluded- volume constant lead to renor- 
malization of correlation length. Substituting Eq.(^) into Eq.(|l2D, one obtains from Eq. (p!0D 
the following expression for the correlation function 

where the renormalized correlation length is determined as 

To obtain Eq.(0) we take the limit g — > in Eq.(0). This approximation is justified for 
consideration of scattering in the optical region because A ^ ^, C,s (see below). Equation 
(0) implies that the renormalized correlation length C,r diverges when the heterogeneity 
parameter g = w'^/v'^ approaches the critical value: 

gc = ^7rc^U^ + Q. (15) 



This means that the system undergoes a second order phase transition at Qc |T3] • Note that 
in O conditions, when f — > 0, the hmit g ^ Qc is reached for any nonzero value of variance 
w ^ because g behaves as Note that for the binary solvent the critical value of $, 
$c, can be found from the equation g{^c) = gc, where g{^) = is determined 

by Eq.(§). 

Now let us present the physical reason for the increases of density fluctuations and the 
renormalized correlation length C,r with the heterogeneity parameter g = w'^/v'^, which 
eventually leads to the phase transition. When is or very small, to maximize the 
entropy (minimize the free energy) polymer chains prefer to take conformations which make 
the density of monomers in solution almost constant with small fluctuations. When g is 
increasing, the solvent is more and more heterogeneous consisting of regions with small and 
large excluded volume constants. In such cases, the monomers prefer the regions with small 
excluded volume constant to minimize the excluded volume energy. This leads to large 
density fluctuations of monomers because of the composition fluctuations of the solvent. 
When the heterogeneity parameter g exceeds the critical value gc, the trend to minimize the 
excluded volume energy having a heterogeneous structure becomes more important than the 
trend to maximize the entropy having almost constant density. Since the polymer solution 
is concentrated, the overlapping between the polymer chains exists and the fluctuations can 
propagate through polymer chains, which leads to the divergence of at gc- 

Now consider how the properties of polymers in a multicomponent solvent just mentioned 
could be reflected on the light scattering. As was stated above the intensity of scattered 
light is proportional to G{q) of Eq. (|T3|) . From Eqs.([T3|), (|1^) and (|15D we obtain 



I(k) (16) 

Here k = ^nsin^/2 with A being the wavelength of the incident light, n being the refractive 
index of the medium, and 6 being the scattering angle. Consider the asymptote of Eq.(|l^). 
Far away from the critical point g <^ gc{^ ^ ^c) the renormalized correlation length is 
almost equal to the bare one ~ C,- In the optical region A ^ ^, neglecting the first term 



in the denominator of Eq. ([T6|) , we find 

/o(fc) ~ e. (17) 

In this case scattered hght is isotropic and the intensity is small because in a good solvent 
the correlation length of a concentrated solution is of the order of a few nanometers 0. In 
the critical region g ~ gc, <^ k except a small region of angles near ^ ^ 0. Therefore for 
scattered light intensity, one has from Eq.(|l6[) 



^(^) ~ 1^ 2 (18) 

loTT'^n^ sm 0/2 

It follows from Eq . (|T8|) that at the critical value g ^ gc of the heterogeneity parameter 
scattered light is strongly anisotropic and the intensity is much larger than that predicted 
by Eq.(p!7D (see below). 

Now we want to apply our theory to polymers dissolved in a binary solvent. The ex- 
perimental system consists of polyethylene oxide (PEO) dissolved in nitroethane+3-methyl- 
pentane (PEO/NE/MP) where nitroethane (NE) is a good solvent for PEO and the methyl- 
penthane (MP) is a poor one. The molecular weight and the concentration of PEO in the 
solvent were 9 x 10^ and 0.075 mg/cc, respectively. Note that a MW = 9 x 10^ molecular 
weight PEO solution in a marginal (mixture of good (NE) and poor(MP)) solvent should be 
treated as a concentrated one because its concentration 0.075 mg/cc significantly exceeds 
the overlapping concentration of the chains (see, for example, [||). The volume fraction 
of MP was controlled in the experiment and the measurements were done at different values 
of the fraction. The refractive indices of NE and MP are very close to each other, 1.39 and 
1.38, respectively. This means that light scattering caused by composition fluctuations of 
the binary solvent could be significant only near the critical temperature of the mixture. 
Note that the binary mixture NE/MP has an upper critical temperature at 26.5° C and the 
critical composition is $o = 65% MP by volume [0 . An argon ion laser with the wavelength 



A = 488nm was used in the experiment. The temperature of the sample used by TKC was 
42° C , far away from the critical temperature. Nevertheless when the volume fraction of 



MP was varied from 0.58 to 0.7 an abrupt increase in the intensity of scattered light was 
observed. The experiment shows more than two orders of magnitude increase in the scat- 
tering intensity within the narrow composition range. At higher composition $ > 0.7 the 
intensity is greatly reduced (see Fig. 1 and Sec. 3 of 0). 

Note that when the fraction of MP is small the average excluded- volume constant, ac- 
cording to Eq.(^, is mainly determined by the good solvent NE, and the heterogeneity 
parameter is small g Qc- Therefore scattered light intensity is determined by Eq.(|l^). 
Increasing the fraction of MP we increase the heterogeneity parameter which reaches the 
critical region g ~ g^. In this case scattered light intensity is determined by Eg . ([T8|) . Com- 
paring the intensities Eqs. ([T7|) and ([T8|) at scattering angle 9 = 90° (as in the experiment), 



we have 

/ / A \^ 



(19) 



Substituting the experimental values |p of A = 488?T,m, n = 1.38 and ^ ~ Anm in Eq. (p!9[), 
we have ///q ~ 100. This is consistent with what was observed in the experiment 0]. Note 
that the other characteristic feature of light scattering: isotropic angular distribution in 
the good solvent regime and anisotropic distribution in the critical region predicted by the 
theory were also observed in the experiment 

Now let us present a consistency check of some experimental data. It follows from Eq. 
(I) that the heterogeneity parameter g = jv^ depends on $ and x = Vp/vg and can be 
written as 

, <l>(l-<l>)(l-x)2 ^ ^ 



Differentiating Eq. (|20| ) with respect to $ and setting the result to 0, we find that the 
maximum of g is reached at $o = l/{x + 1) and gmax = (1 "~ x)'^/Ax. Assuming that 
X ~ 0.4, one has gmax ~ 0.22 and $o ~ 0.7 which is consistent with the critical fraction of 
MP: $ ~ 0.6 — 0.7. Substituting from the experiment the values ^ ~ 4nm, C,s ~ Inm and 
c = 12 X 10^"^ cc~^ into Eq. (|T5|), we find that g^ ~ 0.22, which is consistent with the gmax- 
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Please note that from only one phenomenological parameter x ~ 0.4 we obtain Qmax and $o 
which are consistent with experimental data. 

As mentioned above, the experiment was done at temperatures far above the critical 
temperature of the NE/MP mixture. The question arises of how the scattered light intensity 
will behave when the temperature decreases and approaches the critical temperature of the 
solvent. Obviously the correlation length of the composition fluctuations C,s will increase 
under the temperature decreasing. Accordingly, from Eqs. ([l5|) and (|l^), the scattered light 



intensity will decrease. However, near the critical temperature, although the difference 
between the refractive indices is small, the scattering by the solvent itself becomes essential 
because of the large value of correlation length ^g. Therefore the intensity of scattered 
light will first decrease with decreasing temperature, then begin increasing near the critical 



temperature of the NE/MP mixture. This behavior was observed in later experiments [jl4 

Note that our calculation here corresponds to a mean field approximation. The diver- 
gence of the correlation length at the critical point ~ ^{dc — q) '^ with the exponent 
z/ = 1/2 (see Eq. (|T^) is the reflection of this fact. However in the vicinity of critical point 
g ^ Qc one should take into account large fluctuations. We have shown |]TU[ that effectively 



these fluctuations are described by the field-theoretical 0{N) model at = 0. Therefore 
the correct value of critical index u = 0.588 (see, for example, [l^) follows from this corre- 
spondence. A more detailed experiment on light scattering in the above-mentioned ternary 
system could check this theoretical prediction. 

In summary, we have presented a theory of phase transition due to the heterogeneity 
of the excluded-volume constant in a concentrated polymer solution with multicomponent 
solvent. The abnormal scattering of light observed in the experiment of is a manifesta- 
tion of such critical behavior. Our theory describes the main characteristic features of the 
experiment qualitatively as well as quantitatively. 

We would like to thank Jonathan Dushoff for a critical reading of the manuscript, N. 
Izmailian and Kiwing To for valuable comments. This work was supported in part by the 
National Science Council of the Republic of China (Taiwan) under Contract No. NSC 90- 



2112-M-001-074. 



REFERENCES 

* Electronic address: gevorkia@phys. sinica.edu. tw. 

^ Electronic address: huck@phys. sinica.edu. tw. 

[1] J. des Cloizeaux and G. Jannink, Polymers in Solution, (Clarend Press, Oxford, 1990). 

[2] Y. B. Melnichenko, et al., Phys. Rev. Lett. 79, 5266 (1997); A. Bhattacharya, S. D. 
Mahanti, and A. Chakrabarti, Phys. Rev. Lett. 80, 333 (1998); W. B. Li, et al., Phys. 
Rev. Lett. 81, 5580 (1998); A. M. Hulkarni, et al, Phys. Rev. Lett. 83, 4554 (1999). 
M. A. Anisimov, et al., Phys. Rev. Lett. 85, 2336 (2000). 

[3] H. Fujita, Polymer Solutions, (Elsevier, New York, 1990). 

[4] C. D. Sfatos and E. I. Shakhnovich, Phys. Reports 288, 77,(1997). 

[5] D. Mclntyre and F. Gornick, Light Scattering from Dilute Polymer Solutions, (Gordon 
and Breach Science Publishers, New York, London, 1964). 

[6] P. Debye, Molecular-weight Determination by Light Scattering, p. 123 in 

[7] H. Yawakawa, J. Chem. Phys. 46, 973 (1967) and references therein. 

[8] M. Doi and S. E. Edwards, The Theory of Polymer Dynamics, (Oxford Press, New 
York, 1986). 

[9] K. To, C. A. Kim, H. J. Choi, Physica A 254, 292 (1998). 
[10] Zh. S. Gevorkian and C.-K. Hu, e-print |cond-mat / 020 1 169 . 
[11] R. F. Chang et al, Phys. Rev. Lett. 27, 1706 (1971). 

[12] A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinskii, Methods of Quantum Field 
Theory in Statistical Physics, (Prentice- Hill, Englewood Cliffs, 1969). 

[13] A. M. Wims, D. Mclntyre, F. Hynne, J. Chem. Phys. 50, 616 (1969). 



11 



[14] K. To, private communication, (2001). 

[15] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena{C\axendon Press, Ox- 
ford, 1996). 



12 



